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I. INTRODUCTION 



Most information on the quark and gluon structure of hadrons has come so far from 
inclusive Deep Inelastic Scattering (DIS) type experiments. With an appropriate selection 
of probes and reactions, accurate measurements conducted through the years allowed one 
to map out in detail the different components of proton structure, the Parton Distribution 
Functions (PDFs) in a wide kinematical region of the four- momentum transfer, Q^, and of 
the longitudinal momentum fraction of the proton's momentum, Xbj = Q'^/2Mu, v being 
the energy transfer and M the proton mass. 

Recently, a whole new dimension was added to our understanding of hadronic structure, 
with the observation that in a number of exclusive experiments one can, on one side, study 
a wider range of flavor and spin dependent combinations of PDFs with respect to those ob- 
tained from inclusive scattering, and on the other, new, qualitatively different information 
can in principle be extracted from a specific class of experiments, including Deeply Virtual 
Compton Scattering (DVCS), and hard Exclusive Meson Production (EMP). The informa- 
tion from these processes is coded in terms of "off-forward" contributions, or the Generalized 
Parton Distributions (GPDs) . GPDs allow us to access partonic configurations with 

a given longitudinal momentum fraction, similarly to DIS, but also at a specific (transverse) 
location inside the hadron 0]. 

PDFs are extracted directly from inclusive measurements of the DIS structure functions 
at a given . Perturbative QCD (PQCD) evolution connects the PDF values at with 
the ones at an initial scale, Q\. The initial PDFs are usually given in parametric forms that 
broadly reproduce the behavior expected in a few limiting cases: they include , for instance, 
a Regge-type behavior at low x^j, and a quark-counting type behavior in the limit Xsj 1- 
A number of sum rules such as Adler's, and the momentum sum rule provide additional 
constraints. It is also generally understood that the proton is an "emptier" object dominated 
by its minimal - valence - components at low Q^. The sea quark content is determined both 
by quark- ant iquark pairs and gluons radiations which characterize perturbative evolution 
from the initial low scale, and by "intrinsic" components also expected to be present at 
low scales. PDF parametrizations have become more sophisticated through the years both 
because of the continuous addition of DIS data in increasingly extended kinematical regimes, 
and because of phenomenological developments allowing one to extend the number of hard 
processes from which PDFs can be extracted. It is now possible to describe the proton 
structure functions with relative accuracy in the regime: 10~^ < xbj ^ 0.75, and 1 < < 
10^ GeV^ fl. 

The matching between measured quantities and leading order predictions for DVCS /EMP 
and GPDs should proceed, in principle, similarly to the inclusive case, in view of the factor- 
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ization theorem discussed in Ref. [2j. There are however a few important caveats due to the 
fact that GPDs describe the non-perturbative contribution to an amplitude. Both the real 
and imaginary parts of the amphtude are physical observables in the 7*P —* -yP' process, 
obtained from the interference term for the DVCS and Bethe-Heitler (BH) processes (see 
and reviews in 0, S])- 

The leading order amplitude for DVCS is shown in Fig. [1] along with the relevant kinemat- 
ical variables, namely: the longitudinal momentum fraction taken by the initial quark, X, 

= —qf^, the four- momentum transfer squared between the initial and final proton states, 
t = — A^, and the longitudinal momentum transfer fraction of the initial proton momentum, 
the so-called "skewness" , (. Because of the extra two parameters, ( and t, obtaining initial 
parametrizations in a similar fashion as for inclusive parton distributions is a formidable 
task. Furthermore, experimental measurements of GPDs are remarkably more complicated 
than in inclusive DIS, essentially due to the exclusive nature of DVCS/EMP. Despite the 
high performance of current facilities such as Jefferson Laboratory, one cannot realistically 
expect in the near future a similar amount and quality of data as in inclusive experiments. 

Given the importance of the physics underlying GPDs, and the far reaching consequences 
of their study related to a number of key questions such as the contributions of the partons 
angular momentum to the nucleon's spin, the exploration of 3D spatial images of the nu- 
cleon, and the connection to Transverse Momentum Distributions (TMDs) [7|, [sl], it is most 
important at the present stage, to explore whether currently available inclusive data can 
also provide additional constraints on GPDs. These constraints can both supplement the, 
so far, scarce experimental results obtained directly from DVCS/EMP, and at the same time 
provide a guidance for future precision exclusive experiments at both Jefferson Lab at 12 
GeV, and future colliders. 

Experimental constraints from other-than-DVCS-type data on GPDs are obtained from: 
i) The nucleon form factors providing integrals of GPDs over X at a fixed t; ii) The PDFs, 
representing the 0) limit of GPDs. An additional check is also provided by the 

relation, through simple Fourier transformation, between zero skewedness GPDs and Impact 
Parameter Dependent PDFs (IPPDFs) 0]. Physically meaningful GPDs should in fact 
reproduce the correct behavior of partonic configuration radii obtained from their impact 
parameter space representation. 

Because both the form factors and PDFs are independent of C, these constraints apply 
exclusively to the case = 0, t = — A^. Zero skewness GPDs extracted using available 
experimental constraints were first considered in the initial phenomenological studies of 
Refs. jsl, 0, [nl, [l2]. Fully quantitative fits were subsequently performed in Refs. 0, 
Defining a parametrization at C 7^ 0, however, requires the additional condition of polynomi- 
ality to be satisfied 0] Moreover, a physical interpretation in terms of partonic components 
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becomes less transparent both a.t X = ( - the "stopped returning quark" region - and 
at X < C where the dominating process is scattering from a qq pair emerging from the 
initial nucleon. The only guidance for a parametrization at ( ^ has been provided, so 
far, by the Double Distribution (DD) hypothesis 0, Q, that has a built-in property of 
polynomiality. More recently, the Mellin-Barnes integral representation 17|, and the dual 
representation 18|] were proposed, where the GPDs were obtained within a generalization 
of the anti-Mellin transform approach used for PDFs. Nevertheless, similarly to what found 
for DIS 19, 2o|, the extraction of GPDs from moments formally requires a continuation to 
complex n that, because of an oscillating term inherent in the moments integral, can be a 
source of ambiguities. This has so far hampered an accurate extraction using other than 
simplified models. 

Motivated by this situation, in this paper, we introduce a practical method to extract 
GPDs from experimental data, that can be extended also at C 7^ 0, by using additional 
^-dependent constraints from ab-initio lattice QCD calculations of the first three moments 
of GPDs 2l|, |22|. We stress that differently from the DDs, representing a model calcula- 



tion, our approach is, for the first time to our knowledge, an attempt to obtain a realistic 
parametrization. Given the paucity of current direct experimental measurements of GPDs, 
our goal is to provide more stringent, model independent predictions that will be useful both 
for model builders, in order to understand the dynamics of GPDs, and for the planning of 
future DVCS type experiments. 

Our paper organization is as follows: In Section [III we describe our approach for a physi- 
cally motivated parametrization valid both in the C = and C 7^ cases. We consider the 
unpolarized GPDs, H and E, for which the lattice moments are the most accurate. In Sec- 
tions [Till and IIVI we discuss in detail the C = case: in III we perform a detailed comparison 
with the data on both form factors and PDFs; In IV we show the phenomenology of GPDs 
at C = by illustrating the role of the various quantities: the quarks transverse radii (yq), 
the intrinsic transverse momentum k^, the average X values contributing to the nucleon 
form factors as a function of t. We also discuss the feasibility of the extraction using lattice 
results and introduce our method. For ease of presentation, a number of graphs and more 
quantitative results on the C 7^ case are reported in a following paper 23|]. In Section V 
we draw our conclusions. 



II. A PHYSICALLY MOTIVATED PARAMETRIZATION FOR UNPOLARIZED 
GPDS 

GPDs parameterize the non-perturbative vertex in the DVCS process depicted in Fig. [H 
Scattering from an unpolarized proton (neutron) is described by two independent GPDs: 
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H, and E, from the vector (7^) and tensor {a^y) interactions, respectively, that depend 
on three kinematical invariants, besides the initial photon's virtuality, Q^: the longitudinal 
momentum transfer, C, = Q'^/2{Pq), the four-momentum transfer squared, = —t, and the 
variable X = {kq)/{Pq), representing the Light Cone (LC) momentum fraction carried by 
the struck parton with momentum k. The relations between the variables used in this paper 
and the analogous set of kinematical variables in the "symmetric" system, frequently used 
in the literature are given along with the definitions of the hadronic tensors components in 
Refs. 



2J]. Note the slightly different choice from Ref. j25| for which the sea quarks GPDs 
are also defined at X > (see Appendix A). 

At present, due to both the small experimental coverage mentioned above, and to the 
somewhat less transparent physical interpretation of the observables, an important role in 
the shaping of a parametrization for GPDs is played by the intuition on the underlying 
non-perturbative dynamics. 

In the general case of non-zero skewedness, one distinguishes two kinematical regions: i) 
( < X < 1, where the dominant process is where a quark from the initial proton with LC 
momentum fraction X, is struck by the initial photon and is subsequently reabsorbed in 
the proton (with X — ( > 0); ii) < X < (, where the final quark propagates backwards, 
or, correspondingly, a quark-antiquark pair from the initial proton, with an asymmetric 
partition of LC momenta, participates in the scattering process. 

The physics of the two regions is most easily understood within a field theoretical de- 
scription where the lowest order is given by a covariant quark-nucleon scattering amplitude, 
with the nucleon-quark-diquark vertex being a Dirac matrix multiplied by a scalar function. 
^. Physical intuition on DVCS processes can then be obtained by viewing the covariant dia- 
grams as the sum of all possible time ordered diagrams. The invariant amplitude for DVCS 
corresponds to 4! time ordered diagrams. They are grouped into the following classes of pro- 
cesses (Fig. [2]): (a) all particles moving forward (Fig. [2^); (b) the initial photon splits into 
a quark-antiquark pair that then interacts with the hadronic system; (c) an initial quark- 
antiquark pair originates from the initial proton and scatters from the probe (Fig. Wp)- Each 
process has a corresponding "crossed term" (not drawn in the figure). 



^ The most general case involves a linear combination of Dirac matrices [26l l27l|) 
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We consider the following choice of frame, and four-momentum components:^ 



q = (0;q,0) 
P = (^+7T7T;0,P) 



(la) 
(lb) 



(Ic) 




(Id) 



(Ig) 



(le) 



(If) 



where, moreover, q' = q + A. Eq. (lldp gives the components of a spectator system, i.e. 
a diquark, that is our main assumption in the following sections, namely that the spectral 
distribution of states appearing in principle in the quark correlator can be replaced by 
one state with a given mass. Finally, in the given frame, contributions dominated by the 
hadronic components of the initial photon (case (c)) are absent. More generally, all processes 
where one of the particles is moving backwards (or one of the particles has a longitudinal 
momentum opposite to the protons' one) vanish as inverse powers of P. 

A. The X > ( region 

At X > (, the proton splits into a quark carrying a LC momentum fraction X = /P~^, 
transverse momentum k_|_, and a spectator system with 1 — X,= k^/P~^, and — k_|_ (Fig. [2^). 
Similarly the right side vertex describes the coalescence of the final quark and the spectator 
system into an outgoing proton (all particles are moving forward). With an appropriate 
choice for the P (kkx) and P' — » {k' kx) vertices, i.e. assuming a spectator diquark 
with both scalar and axial vector components, the DVCS matrix element, -F(X, C, ^), can be 

^ We use the notation = (a° = Ea- a^, a^), and = (a° ± a^)/v^. 
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written at leading order in Q^, as (see also Ref. (28|): 

1 



2P^ 



U{P', S') ( i^H\X, C, t) + -^^E\X, C, t) ] U{P, S) 



Vx^/x^ 

l-X 



In Eq. 



1 



VW = — rr{f/(P,5)f/(P',5')7+} 
- = — rr{f/(P,5)I7(P',^')— ^^''A.}, 



(2) 

(3a) 
(3b) 



where the traces over the nucleon spinors are for the same spin, S = S", case. The factors 
Xy/X — are obtained similarly from the traces over the quarks spinors implicit in Eq. ([2]) 



28( 1 . Having taken care of the spin structure, we then model p'^ {k"^, k'"^) as: 

0(A;'2) <p{k^) 



{k\k'^)=M 



_ ^2 ^2 _ ^2 ■ 



(4) 



rUg is the struck quark's mass, 0(/c^) is a scalar vertex function whose form will be specified 
in Section [nil and A/" is a normalization constant. 

By using the components in Eqs. (fTa]) we obtain for H'^{X,(,t): 

d^k^_ A(j){k^,X)<j)*{k'\X) 



H\X,C.t) 



1 - X) [{2Ek){Ep -Ek- Ex)] [{2Ey){Ep, - E^' - Ex)] 

^0(A;2,A)(/>*(A;'2,A) 



MUX) - ^ 



MKX) - j 



-X 



i-x 



(5) 



where we have rendered explicit the connection between the time ordered diagram in Fig. [2^ 
and the covariant expression in Eq. (j4|). Furthermore, in Eq. ([5]): 



Ml{X) = {X- 0/(1 - C)M' -{X- 0/(1 - X)Mf 



and 



A = Af- 



X 



x-C 



(6) 



(7) 



1 -X V X 

The GPD, E''{X, (, t) is modeled similarly to H'^, but imposing a different normalization, 
namely: Me = f^q-Afn, {Mr = TV in Eq. (j4])), where Kg is the quark's q component of the 
anomalous magnetic moment. 



7 



The invariant mass of the spectator, kx = M^, appearing in Eq. ([5]) through Eq. (Q, 
is a flavor dependent parameter. Both in Eq. ([5]), and in the results presented in Section 
mil the value of is considered to be fixed for each configuration. However, a spectral 
distribution in M^^ = {P — k)"^ should in principle be introduced for large values of the 
invariant mass. This affects mainly the low X region, i.e. where is large, and it has 
been successfully reproduced in deep inelastic scattering processes by introducing an 
dependence of the spectral function consistent with Regge behavior [2^. The role of t- 
channel exchanges in DVCS and related processes was addressed recently in 30| where the 
rather extreme point of view was taken that GPDs measure mostly the parton content of 
the reggeons. A full treatment of this important point is beyond the scope of the present 
work and will be considered in a forthcoming paper. Here, we introduce directly a ( and t 
dependent Regge motivated term, in addition to the "diquark" term given by Eq. and 
we study their relative contribution to phenomenology in Section HVl 

In Eq. (|5]) we have written explicitly the dependence on the size parameter A, which is 
similar to the one used in the "overlap representation" based models (cf. e^. the equivalent 
parameter for the more commonly used "gaussian form" discussed in Ref . 3l| ) . We underline 
that the model considered here is, in fact, consistent with the "overlap representation" 
derived for DVCS in 



25 



321, but, due to the covariance of the vertex function, it differs 



from constituent quark models. Due to the covariance of the vertex function, the spectator 
model captures two essential features, or "self- consistency" conditions: i) the GPDs are 
not imposed to be zero at the endpoint X = (, thus allowing for an imaginary part of 
the DVCS amplitude. This is also in accordance with the experimental observation from 
DVCS experiments at both HERMES 33|] and Jefferson Laboratory 3J]; ii) the GPDs are 
continuous at the endpoint X = (. Calculations similar to the one presented here were 
performed both within QED 25|, and in a simplified version of the covariant model with 
scalar particles in 35|. Both cases are however presented as illustrations that are not meant 
to be quantitatively c omp ared to data. On the other side, attempts similar to ours to extract 
GPDs from the data 13| , although physically motivated, are leaning towards mathematical 



forms similar to "PDF-type" parametrizations. 

This paper's goal is to combine both the essential dynamical aspects described above, 
with a fully quantitative analysis that is made possible by the flexibility of our simplified 
model. A direct comparison with inclusive experimental data is only possible in the C = 
region. At ( ^ 0, one needs to include in the analysis the higher moments of GPDs, that 
are (^-dependent, besides the nucleon form factors. Higher moments are currently available 
from lattice QCD 36|] and can be implemented within an extension of our analysis to this 
case. 
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B. A special case: C = 



The case C = where the momentum transfer is entirely transverse, = — A^, plays a 
special role since processes of type (b), and (c) (see Fig. [2]) are suppressed. 
The following relations hold (we set Hi{X,0,t) = Hi{X,t)): 

H'i{X,t = 0)=q{X), (8) 

where q{X) is the parton distribution for quark "g". The transverse DIS structure function, 
Ft{X) = Fi{X) is given by: 

F|(X) = ^H^{X, 0) + ^H\X, 0) + ^H\X, 0) (9a) 

F^iX) = ^H^iX, 0) + ^H\X, 0) + ^H^iX, 0), (9b) 

where we implicitly assume the dependence, and the structure function F2 is obtained 
from Callan-Gross's relation: 2X Fi{X) = F2{X) Furthermore, the following relations: 

[ dXH'^{X,t) = Ff{t) (10a) 
Jo 

[ dXE\X,t) = F^{t), (10b) 
Jo 

define the connection with the quark g's contribution to the Dirac and Pauli form factors. 
The proton and neutron form factors are obtained as: 

^f(2)W = ^^r(2)W - ^Fmit) + ^F^^2)it) (11a) 

Fu2)it) = -^Fr(2)W + IFi%)(^) + \Fm(^), (lib) 



where -Fi(2)(^) "was found to be consistent with zero [37|. In our analysis we fitted linear 
combinations of the integrals of GPDs obtained from Eq. [10] to the electric and magnetic 
form factors, for which the experimental data are more readily accessible: 

Gt\t) = Ff(")(t) + ^F|(")(t) (12a) 
Gl^P\t) = Ff(")(t) + F|(")(t). (12b) 

Eqs. ( Pl9]|10|ll|[T2l) define all the constraints used in our fit. A detailed description of the 
results of the fit is presented in Section IIIII 

C. The X < ( region 

When X < (, the dominating process is the one where a quark with X = k^/P~^, and 
transverse momentum, k_L, and an anti-quark with (—X = k'^ / > 0, and k'_L = A_l — k_L 
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emitted from the initial proton, undergo the electromagnetic interaction (Fig. [2)3). While in 
the calculation of form factors, and of the C = GPDs, process (b) is always suppressed, 
a.t ( ^ 0, it can represent a situation with "all particles moving forward" so long a.s X < (. 
This is evident by inspecting the energy denominators in this kinematical region that are, 
in fact, characterized by similar cancellations as for process (a) aX X > ( (Eq. (Q. 

The physical interpretation of this region still presents, however, a few debatable points. 
Within the overlap representation, GPDs are given exclusively by the higher Fock states - a 
minimum requirement being the (qqqqq) state. The latter are not sufficiently constrained by 



phenomenological studies. Recently, higher Fock states were considered in Ref. [38|] within 
a LC constituent quark model for the pion GPDs, where their contribution was shown to 
be indeed sizable at X = ( (the so-called crossover point). Notice, however, that as X 
decreases, a large number of Fock components would need to be introduced. Quantitative 



calculations were performed both within QED 25|], and in a scalar model 35|, where it was 
shown that the sum of the X > ( and X < ( contributions naturally provides a covariant 
expression, thus satisfying the polynomiality condition. This is the independence of the 
form factor from the parameter (, which if not observed, would signal the presence of an 
artificial frame dependence 32|. An alternative description is where the qq pair contribution 
to GPDs is interpreted as t-channel exchanges 15| (modulo appropriate color factors 32l |) 
of either a single meson, or a tower of mesons, as recently proposed in 



Lacking a uniform picture, and given the important role that will be played by GPDs 
in the X < ( region for the interpretation of a number of experiments: from pp, and pA, 
exclusive reactions 39[, to deep inelastic pion production and other semi-inclusive experi- 
ments at forthcoming facilities, we propose a strategy that can provide further guidance in 
this region. 

Starting from the accurate parametrization of the C = case presented in this paper, 
and obtained directly from experimental data, we subsequently study the constraints for the 
C > case. These are provided, on one side, by the higher order Mellin moments of H and 
E's that govern the behavior with ( of the GPDs. Moments of order n < 3 can be obtained 

22I . On the other side, we notice that the X > region is 



from lattice calculations 



21 



dominated by the same "valence type" configurations as for the C = case described by 
Eq. (Q, and it can be therefore obtained by extending our C = parametrization to this 
kinematics using the same constraints. In this approach, polynomiality is imposed at every 
step, within a "bottom up" type of approach, rather than the "top down" method imphcit 
in both the double distributions [3]. In Section [IV] we illustrate the type of information 
that can be obtained from lattice results in the specific cases of ( = 0. Information on 
the X < ( behavior using higher moments, combined with the experimentally constrained 
X > ( behavior, is extracted according to a deconvolution procedure described in detail in 
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a forthcoming paper 



23| 



III. RESULTS 

We now present our quantitative determination of the unpolarized GPDs, H and E, 
obtained at C = using all available data on the proton and neutron electromagnetic form 
factors, as well as the valence quarks distributions from DIS measurements. Our fit is 



obtained at a low scale, ~ 0.1 GeV^, in line with the approach of Ref. 40| where it is 
assumed that at a low scale the nucleon consists mostly of valence quarks, the bulk of the 
gluon and anti-quark distributions being generated dynamically. 

We reiterate that, although direct measurements of DVCS cannot be currently im- 
plemented, our procedure produces effectively a "parametrization" , in that parameter- 
dependent physically motivated functional forms are fitted to data. The goodness of the 
fit is tested by means of a x^i whose values, along with the parameter errors have been 
quantitatively evaluated and are given below. We, of course, support the future usage of 



DVCS data because they are more directly linked to GPDs 4l| , and we are actively consid 



ering their implementation 23|. The amount of data and their kinematical range is however 



too limited at present to provide sensibly more stringent constraints. We would also like 
to add that parametrization shapes do represent a possible bias in the present analysis as 
well as in any type of fitting (see e.g. discussion in Ref. jl^), for instance the shape of the 
gluon distribution functions has been oscillating through the years between "valence- like" 
to hard-peaked at x — > 0. The problem of the initial bias can be attacked in a similar way 
for GPDs, by tuning in possible new shapes as constraints from new sets of data become 
available, allowing for more refined fitting. 

Starting from Eq. ([5]), with the inclusion of the Regge term discussed in Section III Al 
we obtained two slightly different forms that are both constrained by current experimental 
data: 



Set I 



H\X,t) = G;i;,(X,t)X-"'-^^(i-^)'''* (13) 
E\X, t) = K G^;, (X, t) x^-'-^^'d-^)'^* (14) 



Set II 



H''{x,t) = G;^;;,(x,t)x-""-^"(i-^)''"* (is) 



E''{X,t) = G'|;,(X,t)X-^"-^"(i-^)''"* (16) 
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All parameters except for pi and p2 are flavor dependent; we omit, however, the "q" symbol 
(unless specifically needed) for ease of presentation. The function G has the same form for 
both parametrizations, I and II: 



where 



and 



D(X,k^) = r-m^ (18) 



e = XM^ - -^Ml - (19) 

l-X ^ l-X ^ ' 

= XM' - - (kx-(l-^)A)^ (20) 

l-X ^ l-X ' ^ ' 

m being the struck quark mass, and M, the proton mass. The normalization factor includes 

the nucleon-quark-diquark coupling, and it is set to A/" = 1 GeV^. Eq. (fT7|) . was obtained 
using the following Dirac structure for the vertex in Eq. ([5]): 

T{k, PU = Yl «"(^' ^) ^')' (21) 



where 27 1 



<Pik^X)= 1^3^. (22) 

Finally, the u and d quarks contributions to the anomalous magnetic moments are: 

^^^,-1^'' = 2-03, for g = . . 

" - " ~ \ = 1.67/2, for g = M ' ^^"^^ 

A few comments are in order: 

At present it is important to provide a parametrization that allows one to address a richer 
phenomenology, including the interplay of coordinate and momentum space observables. 
The spectator model used here is ideal because despite its simplicity, it has proven to be 
sufficiently flexible to describe (and predict) the main features of a number of distribution 
and fragmentation functions in the intermediate and large X regions, as well as the uninte- 



grated PDFs 26|, |43| . The spectator system can be a scalar or a spin 1 vector, thus allowing 



^ With this choice of diquark form factor, one ensures that the value of the quark mass, or equivalently the 
position of the pole in the quark propagator, does not play a dynamical role in the model. The values of 
the quark masses are consequently not determined in our fit. 
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us to access both the u and d quark distributions. In Ref. [26|] it was shown that the value 
of the spectator mass in the two cases is crucial in shaping the parametrizations for u and 
d quarks respectively. Here we have let also the mass parameter A be flavor dependent. 

a) Similarly to the case of DIS structure functions the spectator model is not able to repro- 
duce quantitatively the very small X behavior of the GPDs (this proplem is again present in 
the very small t behavior of the nucleon form factors and GPDs). This mismatch is not very 
visible in the results of Refs. 



26 



27 



43| because of the linear scale used in the plots; it is, 
however, responsible for a violation of the baryon number sum rule that becomes particularly 
important in GPD parametrizations since one needs to achieve a precise agreement with the 
nucleon form factors as well. We introduced therefore a "Regge-type" term multiplying the 
spectator model function G^.j^ in E qs. _ ( fT3|14lll5|[T6|) . A similar behavior was considered in 
the "profile functions" of Refs. 
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14| . However, our procedure is distinctively different 
(as discussed also in our results below) because in our case a simultaneous fit to both the 
PDFs from DIS and to the nucleon form factors is performed. In Refs. Q, Q the PDF 
limit (Eq. ([H])) is trivially satisfied, whereas the form factors and the additional constraints 
from the expected Regge behavior are subsequently used to define the GPDs shape. 

in) We considered the two variants shown above, in order to estimate the sensitivity to 
different procedures as also described in ii). Set I and II differ in the determination of E, 
that is in principle unrelated to the forward PDFs, and therefore less constrained by the 
data. 



A. Results of fit from nucleon form factors and PDFs 



The experimental data on the nucleon form factors implemented in the fit are: 



in Ref. 



mG-^lM and G^/G^ 



48 



4^, 



The data selection is the same that was used 



491], where for — t > 1 GeV^ only the measurements based on polarization transfer 



techniques were considered, while the Rosenbluth separation ones were discarded. In the 
fitting procedure all of the form factor data enter simultaneously in the parametrizations 
for H'^ and E"^, respectively. This is at variance with implementing data on the Fi and F2 
form factors, which also require extrapolations from different data sets. By fitting directly 
to the electric and magnetic form factors we obtained a more precise determination since no 
data manipulation is necessary. In particular, an accurate description of the low t region is 
important in view of future comparisons with the lattice determinations 23 1. 

The per number of data points in each data set, as well as for the total number of 
data points, is listed in Table [I] for parametrization Sets I and II. The comparison with 
form factor data is shown in Fig. [3] for the proton, and in Fig. H] for the neutron. The ratio 
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G^/G^^ is shown in Fig. [51 To check to what extent our parametrization is dominated by 
the data from G\,j, we also repeated the fit by weighting each form factor by the number of 
corresponding data points. The effect for the hsted values of "was less than 2%. 



Data Set 


xVNdata Set 1 


xVNdata Set 2 


Data Points 


Gep 


1.049 


0.963 


33 


Gmp 


1.194 


1.220 


75 


GEpjGup 


0.689 


0.569 


20 


Ge^ 


0.808 


1.059 


25 


Gm„ 


2.068 


1.286 


24 


TOTAL 


1.174 


1.085 


177 



TABLE L The x^/^dnta. of the different nucleon form factors obtained from Set I and Set IL Ndata 
is the number of data points available for each set of form factor data. 



The parameters in the t — > limit were determined by fitting to the LO set of Alekhin 
PDFs |50| within the range 10~^ < x < 0.8 and 4 < < 240 GeV^, and imposing the 
baryon number and momentum sum rule such that: 



dXu{X, Q, 



dXd{X, Ql) = 1 

dXX [uiX, Ql) + d{X, Ql 



1. 



(24) 
(25) 
(26) 



The parameters involved in this step, M^, A'' and a'^, q = u,d, obtained at an initial 
scale Ql {Ql = 0.094 GeV^), are listed in Table HTl Notice that: i) they are the same for 
both Sets I and II; ii) in Set I they are by definition the same for the functions H and E (see 
Eqs. (IISIHD)- The parameters in the PDCD evolution were chosen as in CTEQ6L1 5l[ |. 



Flavor 


Mx (GeV) 


A (GeV) 


a 


u 


0.4972 


0.9728 


1.2261 


d 


0.7918 


0.9214 


1.0433 



TABLE IL Parameters fixing the shape of if, q = u,d, at t = 0. The parameters are the same for 
both Set I and Set II. Moreover, they also define the t = limit for E'' in Set I, as it can be seen 
from the definitions given in Eqs. (|13|14p . 



Similar results can be in principle obtained from other current PDF parametrizations 
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52|, however the valence contributions from Ref. 50|] tend to more readily agree with 



the shape given in Eqs. fll3|15l) . The experimental data on DIS structure functions were not 
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used directly, because our parametrization does not include an ansatz for the sea quarks. 
Therefore it is not possible to reproduce within this context the low x behavior of the DIS 
structure function F2{X,Q'^). For the same reason, the estimated errors of Alekhin's PDFs 
were not used in the fit. This aspect is beyond the scope of the present analysis, and will be 
improved in future work on by extending our model to include sea quarks 23|. It should be 
noticed that, as in similar approaches 



26 



27 



35( 1 . higher-order effects might be important 
especially considering the low value of the initial scale, Qq, resulting from the requirement 
that only valence quarks contribute to the momentum sum rule (Eq. fl2^ ). However, on one 



side some of the previous evaluations 27|, |40| do not seem to find quantitatively large higher 



order effects, on the other, the fact that Ql is a parameter in our model, determined itself 



by fitting the valence PDF to the (parametrization of the) data [50| in the large regime, 



lends self-consistency to our procedure (see also [40|, l53|). In order to gauge the effects of 
perturbatuive evolution, a dedicated study of NLO is being considered in upcoming future 
work. 

With the parameters given in Table II the values of the baryon number sum rules for u 
and d are: 1.9998 and 1.0003, respectively. The momentum at the initial scale adds up to 
1.0016. For Set II we imposed the additional normalization condition: 



dXEg{X,t = 0) 



(27) 



the experimental values of k'^, q = u,d being given in Eq. (123|) . while our fitted values were 
1.6715 and -2.0309, for u and d respectively. The comparison with the PDF parametrization 
of Ref. [s^] is shown in Fig. [6l 

The parameters (3i, /52, pi and p2, in Set I, and all parameters defining E in Set II 
(Eq. (|T6|) ). were fitted to the nucleon electric and magnetic form factors, Eqs. (fT2l) . with the 
values of M^, A''', and a"^ fixed as in Table HTl 

In Table IIIII we list the values of the parameters for Set I, with their corresponding 1-a 
errors (which for 6 parameters corresponds to Ax^ = 7.04). The error band on the form 
factors resulting from the 1-a errors on the parameters is displayed in Figs. [3], H] and [5l 



Flavor 


/3i (GeV-2) 


/32 (GeV-2) 


Pi 


P2 


u 


1.9263 ± 0.0439 


3.0792 ± 0.1318 


0.720 ± 0.028 


0.528 ± 0.031 


d 


1.5707 ± 0.0368 


1.4316 ± 0.0440 


0.720 ± 0.028 


0.528 ± 0.031 



TABLE III: Parameters fixing the t behavior of the GPD forms given in Set I (Eqs. (jl3|14p ). The 
subscript 1(2) is for function H'^(E'^), for each flavor. 



In Table HVl we show the values of (32, Pi and p2 for Set II. We do not present their 
corresponding errors because due to the additional number of paramemeters in this variant 
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of the parametrization, the fit tends to be over-determined, and therefore not completely 
quantitative. This problem can be in principle circumvented either with an increased flow 
of new data from DVCS experiments, or by reducing the number of parameters by keeping 
some of their values fixed. While these strategies can be addressed in the future, we consider 
variant II of our fit as an indicative measure of the model dependence of the E component. 



Flavor 


/3i (GeV-2) 


/32 (GeV-2) 


Pi 


P2 


u 


1.9567 


0.1767 


0.742 


0.270 


d 


1.5896 


3.2866 


0.742 


0.270 



TABLE IV: Parameters fixing the t behavior of the GPD forms given in Set II (Eqs. (|15ll6p ). The 
subscript 1(2) is for function H'^(E'^), q = u,d. 

Additional parameters in Eqs. fll5|16l) are given by: = 1.5780 GeV, Mj/^ = 

0.3902 GeV, A^^ = 0.2678 GeV, X^/ = 0.9589 GeV, 5^^ = 0.005381, a^J = 0.7501. Despite 
the larger number of parameters, the overall agreement with the data is not significantly 
different from Set I: the small t shape of the Gm„ data seems is better reproduced, due to 
the increased flexibility, whereas with Set I one obtains a shghtly better description at large 
t. 

B. C = GPDs for u and d quarks 

Our results for the GPDs are presented in Figs. [TIIHl [9l [TOl and [TTl respectively. In Fig. [7] 
we show Hu and plotted vs. X at different values of t for the parameters in Set I (results 
for Set II are within the 2% error band shown in the figure. Plots for and —E^ are 
shown in Fig. [8] for Set I. In Fig. [9] we display results for E^ and E^i, for both Set I and Set 
II. Differently from H, in this case there is a much bigger discrepancy in the shape of the 
curves at the initial scale, Q^, due to the fact that the constraint from the PDFs is missing 
in this case. As is will shown later on, however, PQCD evolution reduces substantially this 
discrepancy. 

In Fig. [To] we show the separate contributions to Hu of the diquark, G\j^, Eq. ( |T71) . and 
Regge terms, R = respectively, at the fixed values of t = —0.08 GeV^, and 

t = —1.8 GeV^. From the figure it appears clearly that the form of the GPDs is determined 
by the diquark shape, with an "envelope" provided by the Regge term. 

We reiterate that we did not attempt at guiding the values for the parameters in our 
Regge motivated terrn, based on results from soft hadron interactions phenomenology, as 
done instead in Ref. [l3| and Ref. [3], respectively. Therefore, our final parametrization, 
does not depend entirely on the Regge or diquark behaviors in either the low or large X 
regions, but on a mixture of both. 
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In Fig. [TT] and Fig. [12] we show the importance of the Regge term, in a quantitative 
fit of PDFs. In Fig. [11] the u-valence distribution, Q^) for our full parametrization, 

including both the Regge and quark-diquark term (see Eqs. (fT3|) . and ([T5]) . is compared to 



both the PDF fit of Ref. 50|], and to fits performed by excluding the Regge term, at = 5 
GeV^. The quality of the "non- Regge" type fits is similar to those performed e.g. in 43|] and 



271] , however it is clear that the diquark fit alone does not provide hard enough distributions 



at low X. This result is also independent from whether one relaxes the constraints from 
the baryon number and momentum sum rules. We reiterate that the term, which is 

essential in obtaining the argument at low X, is also necessary in order to obtain the correct 
value of the baryon sum rule. To test this we performed weighted fits to the low X region 
without the Regge term. By forcing the model to fit accurately at low X resulted in a 
large mismatch at larger X. Similar results are obtained from other variants of the diquark 
model. We conclude that the diquark model is not apt to reproduce the low X behavior, 
and the baryon sum rule which is in turn fundamental for a quantitative fit of GPDs from 
the nucleon form factors. Fig. [12] emphasizes the low X region on a logarithmic scale. It is 
shown in particular, how the final result (full curve) is determined by both contributions of 
our Regge- motivated term (thick dotted curve), and the diquark model term (dashed curve). 
Our fit curve is shown to be in relatively good agreement with the term: where 
a' and /5' take values consistent with Regge determinations of soft hadronic cross sections 
(shaded area in the figure). The interplay between the diquark model and the Regge term 
explains why the values of the "Regge" parameters in our fit differ sensibly from the results 
from hadronic cross sections, nevertheless giving an accurate description of the low X and 
low t regions. 

Both recent parametrizations from Ref. 13] and Ref. [3], are compared with ours for 
j^u{d) [J2] and in Fig. [TU respectively. A similar comparison for is performed in 

Figs. [IS] and [TH] The results presented in the figure were obtained by evolving to Leading 
Order (LO) the parametrizations at the initial scale Qq ~ 0.1 GeV^ to the values of = 1 
GeV^, and = A GeV^, where the parametrizations from 13] and [3] were respectively 
given. We reiterate that in principle higher order effects are important especially considering 
the low values of the scale, and that PQCD evolution for GPDs involves many more subtleties 
that already appear at NLO (see e.g. 3] and references therein). Lacking however, any 
knowledge from experiment, it is still important to take into account the effects of PQCD 
evolution, as indicators of the general trend followed by GPDs. 

At low t we notice a very good agreement in both iJ" and H"^ among all three approaches, 
essentially because the GPDs tend to the forward limit constrained by PDFs, whereas a 
disagreement appears in in the large t behavior. This is clearly an effect of perturbative 
evolution that sets in earlier in our case with respect to Refs. jisl . [l3 ] , and that does not 
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conserve the form of the initial function, as also noticed in [l3l| . One can gain insight on 
this point, as first observed in the preliminary study of Ref. [12], by plotting the quantity: 

^ ^" j;dXH'^iX,t) Fl{t) Jo ^ ' ^' ^ ^ 

representing the average value of X contributing respectively, to the u and d components 
of the nucleon form factors, and Ff (Eq. (HI])). In Fig. Owe plot (^^(t)) at = Q^. 
The contribution to the proton form factor is obtained as: 

Notice that the total momentum carried by the valence u and d quarks is instead given by: 

(X) = 2(X"(0)) + (X'^(O)). (30) 

From Fig. [T71 one can see that at = Q^, {X) = 1, i.e. all the momentum is carried by 
the valence quarks. PQCD evolution implies that: 

{X^it,Q')) = (X^(t,gg)) [as{Q')/aM)Y' , (31) 

where is the strong coupling constant, and the anomalous dimension yields d2 = 0.4267. 
One can therefore observe that although larger values of X tend to dominate the form 
factor as t increases, the effect of PQCD evolution reduces the large X components by a 
Q^-dependent shift. This effect explains the discrepancies in the curves at — t = 5 GeV^ in 
Figs. [T2][T3][T4l[T5l 

A similar behavior is observed for E'^^'^\ although all parametrizations tend to quantita- 
tively differ also at low t, since they are not constrained by the PDFs. 



IV. PHENOMENOLOGY 

With the results of a precise fit at hand, we can address some issues in the phenomenology 
of GPDs at zero skewness. 

A. Coordinate Space Observables 

A most interesting aspect is the relation with the nucleon Impact Parameter dependent 
PDFs (IPPDFs), q{X,h), related via Fourier transformation to H'^{X,t = -A\) [J]: 

q{X,h) = j d^Ae^^-^H^X^t), (32) 
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q{X, b) is the probability of finding a quark in the proton carrying momentum fraction X, 
at impact parameter b. 

The quark's average impact parameter can be derived as: 



j(Ph g(X,b) 



A-\ogH'^{X,t) 



(33) 



from which an average "interparton distance" 5J] can be defined as: 



{yl{x)) 



mx)) 



(i-x)- 

Similarly, E'^ can be interpreted in terms of a distribution function in a transversely 
polarized target, through: 

g^(X,b) = g(X,b) - d'Ae^'^-''E^{X,t), (35) 

where polarization is along the x axis and by is the component of b along the y axis. 
q-^{X,h) measures the probability of finding a quark carrying momentum fraction X in 
the transversely polarized proton, at impact parameter b. The average shift in the quark's 
distance along the y-axis for polarization along the x-axis direction is obtained from Eq. fl35|) 

Jd^h q^{X,h)by _ 1 E«(X,0) 



as 



mx)) 



(36) 



(37) 



Jd^hq^{X,h) 2MHi{X,0)' 
The shift relative to the spectator quarks is obtained analogously to Eq. ( !34l) as: 

{sq{X)) = YZTjf- 

The average interparton distances, and the transverse shifts are shown in Fig. [18], and Fig 
respectively. 

In our approach the radii are a result of the fitting procedure rather than an additional 
constraint as in Refs. [isl, Q- 

We studied the role of the Regge-type and diquark terms 
from Eq. ([5]). Because of the factorized form, {y'^{X)) can be expressed in fact as the sum of 
the two terms. We find a larger than intuitively expected contribution of the diquark term 
to both the d and u quarks interparton distances. Notice that the interparton distances are 
subject to PQCD evolution. In Fig. [18] we display results at our initial low scale, Q^. The 
effect of evolution is shown by plotting the total interparton distance - including both Regge 
and di quarks terms - at = 4 GeV^. One can see that the interparton distances tend to 
decrease with Q^, as the term dH/dt in Eq. (!33|) evolves more steeply than H. 
The total radius squared of the nucleon is obtained by considering: 



in ) 



j dX q{X){bl{X)) 



t=0 



4|logFf(t) 



(38) 



i=0 
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and using isospin symmetry: 

(rf) = l{rf)-l{rr) (39a) 

{rf) = + (39b) 

The charge radii are given by: 

(-'e) = (rf) + (40a) 

{<) = (rf) + (40b) 

We find: (r]'^) = 0.654 fm^, {rf) = 0.666 fm^, {rf) = 0.650 fm^ (r^^^ = 0.0078 fm^. The 
calculated values for the total charge radii are: (r^) = 0.76 fm^ and (r^^) = —0.118 fm^, 
in agreement with the experimental results. Therefore our results for the fit parameters 
«i 1 — 1.2 in Eq. ([5]) are in line with the constraints on the Regge parameters studied in 
Ref. Q. 

The transverse shift, Sq, (Fig. [T9l) constitutes in principle a test on the GPD E"^. Our 
results, plotted at the initial scale, show a marked difference between Set I and Set II. One 
should keep in mind however, that PQCD evolution largely diminishes the discrepancy. 



B. Intrinsic Transverse Momentum 



An observable related to transversity that can be accessed witihin our model, is the 
partons' average intrinsic transverse momentum, {k'j_{X))q. Notice that k_|_, although not 
Fourier conjugate to b, it can be related to q{X, h) as shown in Ref. 12|. In other words, 
k_|_ is not directly observable in DVCS type processes. However, it can be evaluated using 
the same input to Eq. as: 



(4(A-)>, = '\ (41a) 





0(P,A2)| 






1 


2 



{kliX)) ^ {kl{X))p = ^{kliX))^ + likliX)), (41b) 

In Fig.[20]we show {k'j_{X)) from our model. In order to assess the range of possible variations 
or this observable, we compare our evaluation with the values extracted from Refs. [3] and 
1^ , by assuming a gaussian k± dependence of the vertex functions that could originate such 
parametrizations. We also compare with the hypothesis originally advanced by Burkardt 
|55| on the form of the combined X and t dependences of the gaussian's exponent (or the 



Such relation is of general validity, unrelated to non-relativistic many body theory as quoted in 
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profile function in [l3|, Il4| ' 
(SIDIS) parametrizations 
a few guide 
reactions 



. Finally, we compare with the values used in Semi-Inclusive DIS 



561] . The exploratory work presented here is aimed at defining 



ines for future quantitative studies of the connection between GPDs and SIDIS 



STj . In Fig. [21] we show the ratios: 



Rkl(t) = Fi(2)(t,A;„ax)/Fi(2)(t) (42) 
RflM = i^i(2)(t,X„,,)/Fi(2)(t) (43) 

in order to determine what k±, and X components the form factors are dominated by. The 
numerators were obtained by setting the upper limit of integration in Eq. ( fTOl) to different 
values of k± = kmax, and X = Xmax, respectively. It can be clearly seen that for all values 
of t the form factor ratio, -R[,2^^(t), is saturated by setting kmax ~ 1 GeV (for F2, even 
larger values of k± seem to be important). This is in turn an indication of a semi-hard 

distribution, as opposed to the soft gaussian forms used elsewhere. On the other side, the 

1(2) 

ratio RxJ^^{t) clearly shows the coupling between the X —> 1 behavior of the GPDs and 
the large t behavior of the form factors. 

C. Implementation of Lattice QCD Results 

Lattice QCD provides the only "model independent" constraints that are necessary to 
parametrize GPDs at C 7^ 0. It is, however, important to illustrate the type of information 
that can be obtained on GPDs starting from the C = case. The C 7^ case requires in 
addition, a more involved "deconvolution" procedure from the first three moments that will 



be described in a following dedicated paper [23 1 



Current lattice results can reproduce the dipole fall-off of the form factors up to — t ~ 3 
GeV^. However, even after performing a linear extrapolation to low values of the pion mass, 



such calculations overshoot the experimental results [13|, [21[, [22]. As a result, predictions 
from the lattice are at the moment characterized by a rather large uncertainty. This be- 
comes problematic especially for the higher moments of GPDs where no comparison with 
experimental results can be made. In order to evaluate the impact of such uncertainty on 



possible extractions of GPDs, we performed a test based on a prescription proposed in 36] 
by Schierholz, according to which the value of the dipole mass, A„, appearing in: ^ 

Mn{t)= / dXH{X,t)X^-' = Mn{0)- (44) 

Jo '^+^^ 



Wc reiterate that Eq. (|^^ is not of general validity but a result of fits to the lattice calculations of 2l|, . 
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can be first extracted from tlie lattice evaluations for n < 3, and subsequently extended to 
all n values by performing a fit based on a Regge motivated ansatz of the type: 



(in Eq. dBD iJ = W - H'^ = H"-'^ evaluated at = 4 GeV^). The n-dependence of the 
moments described above allows one to perform the anti-Mellin transform, thus obtaining 
GPDs, analytically. Results using Eqs. fl44|45p . are compared with the evaluations from our 
analysis aX = A GeV^, in Fig. [221 The band in the figure includes the estimate of the 
lattice error. One can observe a good agreement with both our results and other current 
parametrizations for the at t < 2 GeV^. The agreement however, deteriorates at lower t 
for the d-quark. We therefore conclude that within the range of t shown in the figure, it is 
acceptable to use current lattice results in order to extend our analysis of the extraction of 
GPDs at C 7^ 0. 

V. CONCLUSIONS AND OUTLOOK 

In this paper we proposed a fully quantitative physically motivated parametrization of 
generalized parton distributions that is constructed directly from a covariant model for the 
quark-nucleon scattering amplitude. Therefore, we do not implement a specific form for the 
forward limit given by the parton distribution functions of DIS, but we obtain that limit 
from fitting directly to the valence contribution to DIS data. This allows us in particular 
to better study the role of Regge-type exchanges, that are disengaged, in our case, from the 
specific form of parton distributions. The other constraints defining our parametrization at 
zero skewness are provided by the electric and nucleon form factor data. 

The advantages of this approach are that on one side, using the same initial formalism, 
we can predict additional quantities such as the the unintegrated, /c^-dependent, parton 
distribution functions. This degree of flexibility is desirable in view of both future interpre- 
tations of both coordinate space observables, and of transversity. Furthermore, analyzing 
directly the vertex structure of the scattering amplitude allows us to easily extended our 
predictions to the non-zero skewness case. This can be, in fact, obtained by extending the 
range of the kinematical variables in our expressions for the quark scattering dominated 
region, and by considering separately the X < ( region where a quark-antiquark pair from 
the initial proton participates in the scattering process. Additional constraints need however 
to be provided by the higher, (^-dependent, moments of the GPDs. These are in principle 
available from recent lattice calculations but are necessarily fraught with uncertainties. Here 
were able to assess the impact of such ambiguities on the extraction of GPDs. We conclude 
that current extrapolated lattice values can be used for t < 2 GeV^. A detailed description 




(45) 
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of the ( dependent parametrization is given in a forthcoming paper. 

With this type of parametrization in hand we can on one side provide predictions for 
both recent 



41 



58 



and future DVCS measurements at Jefferson Lab. On the other, our 
approach is geared towards providing a practical and more flexible method to reconstruct 
generalized parton distributions from their first few moments. 
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APPENDIX A: CONNECTION BETWEEN DIFFERENT SETS OF KINEMAT- 
ICAL VARIABLES 



By defining [2|: x = + k'+)/{P+ + P'+), ^ = A+/(P+ + P'+), and t, the following 
mappings with the kinematical variables defined in the paper obtain: The regions < x < 
^ and ^ < X < 1 are mapped into < X < ( and ^ < X < 1 respectively, via: 



X 

C 



while the region — 1 < x < — ^ maps into ( < X < 1, via : 

-x + ^ 



X 

C 



(Ala) 
(Alb) 

(A2a) 
(A2b) 



APPENDIX B: ANALYTIC EXPRESSIONS FOR THE NUCLEON'S RADIUS 



We give the analytic expressions for (6^(X)), obtained from the factorized form in 



Eq. ([13 



d 



m^)) = 4-iogi/''(x,t) 



i=0 



1 

2A 



dA ' dA 

mX))G+{bl{X))n) 
3 1 -X 



A=0 



5 (-M2 + - a: 



-(3,{i-xy^ hgx 



with A 



t, Gl, given in Eq.([I7D, and P = x-^-P^'-^-^^^'K 



(Bl) 



[1] D. Muller, D. Robaschik, B. Geyer, F. M. Dittes and J. Horejsi, Fortsch. Phys. 42, 101 (1994) 

[2] X. D. Ji, Phys. Rev. D 55, 7114 (1997) 

[3] A. V. Radyushkin, Phys. Rev. D 56, 5524 (1997) 

[4] M. Burkardt, Int. J. Mod. Phys. A 18, 173 (2003); ihid Phys. Rev. D 62, 071503 (2000) 

[Erratum-ibid. D 66, 119903 (2002)]. 
[5] W.-M. Yao et al., J. Phys. G 33, 1 (2006), "2006 Review of Particle Physics". 
[6] A. V. Behtsky D. Mueller and A. Kirchner, Nucl. Phys. B 629, 323 (2002) 



24 



[7] M. Diehl, Phys. Kept. 388, 41 (2003). 

[8] A. V. Belitsky and A. V. Radyushkin, Phys. Kept. 418, 1 (2005) 

[9] P. Stoler, Phys. Rev. D 65, 053013 (2002) 

[10] A. Gardestig, A. P. Szczepaniak and J. T. Londergan, Phys. Rev. D 68, 034005 (2003) 

[11] A.V. Afanasev, |hep-p h/980 8291 t [hep"-ph/9910565 [ 



[12] S. Liuti and S. K. Taneja, Phys. Rev. D 70, 074019 (2004) 
[13] M. Diehl, T. Feldmann, R. Jakob and P. Kroll, Eur. Phys. J. C 39, 1 (2005) 
[14] M. Guidal, M. V. Polyakov, A. V. Radyushkin and M. Vanderhaeghen, Phys. Rev. D 72, 
054013 (2005) 



[15] A. V. Radyushkin, arXiv:hep-ph/0101225 



[16] M. Vanderhaeghen, M. Guidal and P. Guichon, "VGG code", private communication. 

[17] D. Mueller and A. Schafer, Nucl. Phys. B 739, 1 (2006) 

[18] V. Guzey and M. V. Polyakov, Eur. Phys. J. C 46, 151 (2006) 

[19] G. C. Fox, Nucl. Phys. B 131, 107 (1977). 

[20] F. J. Yndurain, Phys. Lett. B 74, 68 (1978). 

[21] M. Gockeler et al. [QCDSF Collaboration], Phys. Rev. Lett. 92, 042002 (2004); M. Gockeler 
et al, Nucl. Phys. Proc. Suppl. 153, 146 (2006) [arXiv:hep-l at/0512011] 

[22] P. Hagler, J. Negele, D. B. Renner, W. Schroers, T. Lippert and K. Schilling [LHPC collabo- 
ration], Phys. Rev. D 68, 034505 (2003). 

[23] S.Ahmad, H.Honkanen, S. Liuti and S.K. Taneja, "Generalized Parton Distributions from 
Hadronic Observables: ^ 7^ 0, in preparation. 

[24] K. J. Golec-Biernat and A. D. Martin, Phys. Rev. D 59, 014029 (1999) 

[25] S. J. Brodsky, M. Diehl and D. S. Hwang, Nucl. Phys. B 596, 99 (2001) 

[26] H. Meyer and P. J. Mulders, Nucl. Phys. A 528, 589 (1991). 

[27] W. Melnitchouk, A. W. Schreiber and A. W. Thomas, Phys. Rev. D 49, 1183 (1994) 

[28] S. Liuti and S. K. Taneja, Phys. Rev. C 72, 034902 (2005) 

[29] S. J. Brodsky, F. E. Close and J. F. Gunion, Phys. Rev. D 8, 3678 (1973). 

[30] A. P. Szczepaniak and J. T. Londergan, 'arXiv:h ep-ph/0604266, 

[31] J. Bolz and P. Kroll, Z. Phys. A 356, 327 (1996) 

[32] M. Diehl, T. Feldmann, R. Jakob and P. KroU, Nucl. Phys. B 596, 33 (2001) [Erratum-ibid. 
B 605, 647 (2001)] 

[33] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 87, 182001 (2001) 
[34] S. Stepanyan et al. [CLAS Collaboration], Phys. Rev. Lett. 87, 182002 (2001) 
[35] S. J. Brodsky and F. J. Llanes-Estrada, Eur. Phys. J. C 46, 751 (2006) 

[36] G.Schierholz, in GPD 2006: Workshop on Generalized Parton Distributions, June 2006, ECT* 



Trento, Italy - http://gpd.gla.ac.uk/gpd2006/index.php; and private communication. 



[37] A. Acha et al. [HAPPEX collaboration], ar Xiv:nucl-e x/0609002 . 

[38] C. R. Ji, Y. Mishchenko and A. Radyushkin, Phys. Rev. D 73, 114013 (2006); H. M. Choi, 
C. R. Ji and L. S. Kisslinger, Phys. Rev. D 64, 093006 (2001) 



[39] V. Barone, et al, hep-ex 0505054 



[40] M. Gliick et al, Z. Phys. C 41 (1989) 667; ibid Z. Phys. C 48 (1990) 471. 



[41] C. Munoz Camacho et al. [Jefferson Lab Hall A Collaboration], arXiv:nucl-ex/0607029 
[42] J. Pumplin, AIP Conf. Proc. 792, 50 (2005) 

[43] R. Jakob, P. J. Mulders and J. Rodrigues, Nucl. Phys. A 626, 937 (1997) 
[44] L. E. Price, J. R. Dunning, M. Goitein, K. Hanson, T. Kirk and R. Wilson, Phys. Rev. D 4 
(1971) 45; G. G. Simon, C. Schmitt, F. Borkowski and V. H. Walther, Nucl. Phys. A 333, 



25 



381 (1980). 

[45] G. Hohler, E. Pietarinen, I. Sabba Stefanescu, F. Borkowski, G. G. Simon, V. H. Walther and 
R. D. Wendling, Nucl. Phys. B 114, 505 (1976); E. J. Brash, A. Kozlov, S. Li and G. M. Huber, 
Phys. Rev. C 65, 051001 (2002) 

[46] T. Eden et ai, Phys. Rev. C 50, 1749 (1994); I. Passchier et al, Phys. Rev. Lett. 82, 4988 
(1999); C. Herberg et al, Eur. Phys. J. A 5, 131 (1999); J. Golak, G. Ziemer, H. Kamada, 
H. Witala and W. Gloeckle, Phys. Rev. C 63, 034006 (2001); R. Madey et al. [E93-038 
Collaboration], Phys. Rev. Lett. 91, 122002 (2003); G. Warren et al. [Jefferson Lab E93-026 
Collaboration], Phys. Rev. Lett. 92, 042301 (2004); R. Schiavilla and I. Sick, Phys. Rev. C 
64, 041002 (2001); D. Rohe et al, Phys. Rev. Lett. 83, 4257 (1999). 

[47] G. Kubon et al, Phys. Lett. B 524, 26 (2002); H. Anklin et al, Phys. Lett. B 336, 313 
(1994); H. Anklin et al, Phys. Lett. B 428, 248 (1998); A. Lung et al, Phys. Rev. Lett. 
70, 718 (1993); W. Xu et al, Phys. Rev. Lett. 85, 2900 (2000); W. Xu et al [Jefferson Lab 
E95-001 Collaboration], Phys. Rev. C 67, 012201 (2003); S. Rock et al, Phys. Rev. Lett. 49, 
1139 (1982). 

[48] M. K. Jones et al [Jefferson Lab HaU A Collaboration], Phys. Rev. Lett. 84, 1398 (2000); 
O. Gayou et al [Jefferson Lab HaU A Collaboration], Phys. Rev. Lett. 88, 092301 (2002); 
T. Pospischil et al [Al Collaboration], Eur. Phys. J. A 12, 125 (2001); B. D. Milbrath et al 
[Bates FPP collaboration], Phys. Rev. Lett. 80, 452 (1998) [Erratum-ibid. 82, 2221 (1999)] 

[49] J. J. Kelly, Phys. Rev. C 70, 068202 (2004). 

[50] S. Alekhin, Phys. Rev. D 68, 014002 (2003) 

[51] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207, 
012 (2002) 



[52] R. S. Thorne, A. D. Martin and W. J. Stirling, arXiv:hep-ph/0606244 



[53] A. W. Schreiber, A. W. Thomas and J. T. Londergan, Phys. Rev. D 42, 2226 (1990). 

[54] D. E. Soper, Phys. Rev. D 15, 1141 (1977). 

[55] M. Burkardt, Phys. Rev. D 66, 114005 (2002). 

[56] M. Anselmino, M. Boglione, U. D'Alesio, A. Kotzinian, F. Murgia and A. Prokudin, Phys. 

Rev. D 72, 094007 (2005) [Erratum-ibid. D 72, 099903 (2005)] 

[57] M. Diehl and Ph. Hagler, Eur. Phys. J. C 44, 87 (2005). 

[58] H. Avakian, private communication. 



26 




FIG. 1: Amplitude for DVCS at leading order in . The light-cone coordinates for the active 
quarks and nucleons are explicitly written. 




FIG. 2: Time ordered diagrams for DVCS: (a) dominant contribution in X > region; (b) a qq 
pair is first produced from the nucleon and subsequently interacts with the photons. This process 
dominates the X < region; (c) the initial photon splits into a qq pair that interacts with the 
hadronic system. The crossed-terms where two of the particles in the same class are switched, are 
not shown in the figure. 
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FIG. 3: The proton magnetic and electric form factors, G^r and Gp, respectively, divided by the 



dipole, Gd = 1/(1 + (5V0-71GeV2)2, plotted vs. A^. Experimental data from [4J1 {Ge^)\ [451 
{Gmp )■ The full line was obtained using parametrization I, Eqs. ()13ll4p : the dot-dashed line 
corresponds to parametrization II, Eqs. (I15|16p . 




[GeV^] 



FIG. 4: The neutron magnetic form factor, GJj, divided by the dipole, Gd = 1/(1 + 
A^/0.71 GeV^)^, and electric form factor, G^;, respectively, plotted vs. A^. Experimental data 
from 



i7|](GM„) and [40] {G%). Notation as in Fig. [3l 
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FIG. 5: The ratio of proton electric and magnetic form factors, HpG^j{t)/G^{t). Experimental 
data from [43] . Notation as in Fig. [3l 
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FIG. 6: (color online) Parton distribution functions, Xuy{X) and Xdj]{X) plotted vs. X, at the 
initial scale, Qq = 0.3 GeV, and at Q = 5 GeV. The parton distributions parameters, including 
the value of the initial scale were obtained directly from our fit. The fitted parameters are: Mx, 
A, a for both parametrization I and II. Our results are shown together with the LO set of Alekhin 
PDFs 



50(] used in the fit. 
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FIG. 7: The generalized parton distributions, (left panel) and (right panel) obtained from 
Parametrization I, defined by Eq. (fT3|) . plotted vs. X at —t = 0,0.1,0.3,1,10 GeV^, respectively, 
at the initial scale Qq = 0.09 GeV^. Results for Parametrization II, Eq. (jlSp are very similar. They 
are included within the 2% band displayed in the figure 
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FIG. 8: The generalized parton distributions, (left panel) and —E^ (right panel) obtained from 
Eq. (|14p - Parametrization I - plotted vs. X aX, —t = 0,0.1,0.3,1,10 GeV^, respectively, at the 
initial scale = 0.09 GeV^. 
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FIG. 9: Left panel: obtained with Parametrization I (full lines) and II (dashed lines), respec- 
tively, for different values of -t: -t = 0,0.1,0.3,1 GeV^, at the initial scale Ql = 0.09 GeV^. 
Right panel: the same for —Efi . 
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FIG. 10: Contributions to Hu of the diquark term: G^^^, Eq. (fT7|) . and of the Regge terms: 
R = ^-a-/3(i-X)fu^ respectively, at t = -0.08 GeV^ (top) and t = -1.8 GeV^ (bottom). From 
the figure one can see that the form of the GPDs is determined by the diquark shape, with an 
"envelope" provided by the Regge term. 
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FIG. 11: (color online) Role of the Regge-motivated term, R = X " f^^^ x)nt^ jj^^ order to 
accomplish a quantitative fit of the PDF Uv{X,Q'^). 
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FIG. 12: (color online) Role of Regge-motivated term in Eqs. (113|15p . Thick dotted line: 
^-a-p{i-x)pit^ with fit parameters from Tables II and III; dashed line: quark-diquark term, ; 
Full line: total value, Eqs. p3ll5p . Shaded area: Regge term: with fit parameters from 

phenomenological Regge fits. All curves were calculated at t = —0.1 GeV^. 
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FIG. 13: (color online) Comparison with the quantitative extraction of GPDs from data from 
Ref. [l^. Hu (left) and Hd (right) as a function of X for —t = 0, 0.3, 5 GeV^, for our Parametriza- 
tions I and II, respectively evolved at = 4 GeV^ used in Ref. IJ] While a clear agreement is 
seen at low values of t, this becomes worse at larger values of t. 
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FIG. 15: (color online) Comparison with the quantitative extraction of GPDs from data from 
Ref. Eu (left) and —Ed (right) as a function of X for —t = 0, 0^, 5 GeV^, for our Parametriza- 



tions I and II, respectively evolved at = 4 GeV^ used in Ref. 13|. 
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FIG. 17: (color online) Average value of X, Eq. (j28l plotted vs. —t at = Qq for the u and 
d quarks contributions to the proton Dirac form factor, Eq. (jlip . The average vale of X for the 
proton is also shown at = 4 GeV^ (dot-dashed line). 




FIG. 18: (color onhne) {yl{X)) (left panel) and {yj{X)) (right panel) in our model, plotted vs. 
X. The contribution of the Regge and diquark term, respectively (see text) are shown separately 
along with the total result, at the initial scale. The total result is then evolved to Q'^ = 4 GeV^. 



38 



A 1-2 

R 

r= 1 

V 

0.8 
0.6 



0.4 



0.2 



Set I 
Set II 




','T;vt-,j...,....j...,., , 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



1.2 



0.8 



0.6 



0.4 



0.2 



— Set I 




Set II 











0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 19: (color online) {su{X)) (left panel) and {S(i{X)) (right panel) in our model, plotted vs. X. 
All curves are at the initial scale, Qq. The average interparton distance: (?/^(X))^/^ is shown for 
comparison. 
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FIG. 20: (color online) The average intrinsic transverse momentum, Eq. ()41bp . compared to the 
value extracted from GPD parametrizations (Refs. [l^, 14]); to the t-dependence conjecture of 



Ref. 



55( 1 ■ and to models used in Semi- Inclusive DIS (SIDIS), Refs. [43| and [561 ]. 
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FIG. 21: Contribution of both the intrinsic transverse momentum components (upper panels) and 
X components (lower panels) to the proton form factors, Ff{t) (left), and -Ff (0 (I'ight). 
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FIG. 22: (color online) Comparison of both the results of our analysis, and of Ref. 13], with an 
extraction of GPDs from lattice calculations according to the prescription of Ref. 36|]. The band 
includes an estimate of the error from lattice calculations. 
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